Search for the Lepton Number Violating Process 
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The NuTeV experiment at Fermilab has used a sign-selected neutrino beam to perform a search 
for the lepton number violating process z7 M e~ — > /i~F e , and to measure the cross-section of the 
Standard Model inverse muon decay process z^ M e~ — » fi~v e . NuTeV measures the inverse muon 
decay asymptotic cross-section o/E to be (13.8 ± 1.2 ± 1.4) x 10 -42 cm 2 /GeV. The experiment also 
observes no evidence for lepton number violation and places one of the most restrictive limits on 
the cross-section ratio cr(F M e _ — » /i~V e )/<j(u^e~ — > ^,~v e ) < 1.7% at 90% C.L. for V— A couplings 
and < 0.6% for scalar couplings. 



Neutrino-lepton interactions provide an excellent tool 
to study the properties of the weak interaction. Such 
purely leptonic processes experience no interference from 
strong coupling terms, and thus provide a direct channel 
to investigate the nature of the weak force. The inverse 
muon decay (IMD) process: 
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allows one to make an accurate determination of the 
vector/ axial- vector (V— A) nature of the weak interac- 
tion Q. This process is also sensitive to scalar couplings 
and right-handed currents. 

An experiment with separate neutrino and anti- 
neutrino beams can search for the process: 
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Such an interaction is forbidden by the Standard Model, 
since it violates lepton family number conservation 
(AL e = — ALp = 2). Theories which incorporate multi- 
plicative lepton number conservation H, left-right sym- 
metry H , or the existence of bileptons [Q] allow for such 
processes to occur. 

The NuTeV neutrino experiment at Fermilab has 
investigated these processes in its high-energy, sign- 
selected neutrino beamline. Although the NuTeV in- 
verse muon decay measurement is dominated by system- 
atic uncertainties, the search for lepton number violation 



FIG. 1. Schematic of the NuTeV detector, showing the 
calorimeter and the toroid spectrometer. 



(LNV) processes is very sensitive because the relevant 
backgrounds are highly suppressed. 

The experiment collected data during the 1996-1997 
fixed target run, receiving a total of 2.9 x 10 18 800 GeV 
protons striking a BeO target. Pions and kaons pro- 
duced in the interaction were focused using the Sign- 
Selected Quadrupole Train (SSQT) § and aimed toward 
the NuTeV detector at a 7.8 mrad angle relative to the 
primary proton beam direction. The SSQT enabled the 
detector to be exposed to either pure neutrino or pure 
anti-neutrino beams. NuTeV received 1.3 x 10 18 and 
1.6x 10 18 protons on target for neutrino and anti-neutrino 
running modes, respectively. The fractional contamina- 
tion from wrong-sign meson decays was below 5 x 10~ 3 
||. Pions and kaons decay to neutrinos as they travel 
through a 440 m vacuum pipe; undecayed hadrons are 
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filtered out in a beam dump at the end of the pipe. The 
neutrinos pass through about 900 m of earth berm shield- 
ing before reaching the NuTeV neutrino detector. 

The NuTeV detector Q, located 1.4 km downstream 
of the primary target, consists of a segmented iron- 
scintillator sampling calorimeter, followed by a toroid 
spectrometer (see Fig. |l|). The calorimeter is composed 
of 42 segments, each segment consisting of four 2 inch 
thick steel plates, two liquid scintillator counters, and 
one drift chamber. The calorimeter serves as a neutrino 
target with a fiducial mass of 350 tons. The scintilla- 
tion counters measure the deposited hadronic energy and 
the drift chamber determine the position and direction of 
the outgoing muon. The toroid spectrometer uses a 15 
kG toroid magnetic field to measure the charge and en- 
ergy of muons exiting from the calorimeter. The toroid 
magnetic field is configured so as to always focus muons 
coming from the selected neutrino beam for neutri- 
nos, [i + for anti- neutrinos). The energy resolution and 
response of the detector is measured directly using a sep- 
arate beam of hadrons, muons, and electrons at varying 
energies. The hadronic energy resolution of the calorime- 
ter is a jE = (0.024 ± 0.001) © (0.874 ± 0.003)/V£, 
and the electromagnetic energy resolution is a/E = 
(0.04 ± 0.001) © (0.52 ± 0.01)/ VE @. The resolution of 
the muon energy as determined by the toroid spectrom- 
eter is Ap/p = 11%, limited predominantly by multiple 
scattering. 

The selection criteria for the inverse muon decay mea- 
surement and the lepton number violation search were 
similar, since the characteristic signatures of the pro- 
cesses are nearly identical. Candidate events were se- 
lected based on the following criteria: the event occurred 
during the beam gate, had its interaction vertex within 
the fiducial volume, and had a single \x~ reaching the 
toroid spectrometer. The muon was required to be well 
contained within the toroid and to have an energy be- 
tween 15 and 600 GeV. The muon angle was also re- 
quired to be less than 150 mrad with respect to the beam 
axis. To reduce the number of cosmic ray muons enter- 
ing the selection sample, events which contained signifi- 
cant activity upstream of the reconstructed vertex were 
removed. The hadronic energy of the interaction was 
required to be less than 3 GeV. Finally, the neutrino 
beam running mode determined the sample into which 
the events were placed. For IMD, we required a /i _ 
in neutrino mode; for LNV candidates, a fj,~ in anti- 
neutrino mode. For the LNV sample, we also placed 
an additional requirement on the x 2 of the muon track 
within the toroid, in order to minimize events where the 
charge of the muon was misidentified. 

Because both IMD and LNV events involve neutrino 
scattering from an electron, there exists a kinematic limit 
on the transverse momentum of the muon: pi < 2m e £ , „, 
where m e is the mass of the electron and E v is the neu- 
trino energy. We therefore apply an energy-dependent 
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FIG. 2. Transverse momentum distributions for data 
(crosses) and Monte Carlo (solid) for right-sign neutrino 
events (left) and right-sign anti- neutrino events (right). The 
plot on the left is broken down into background only (dashed) 
and IMD signal (dotted). The plot on the right shows Monte 
Carlo background only. 



requirement on the transverse momentum of the event in 
order to further isolate signal events. The cut requires 

Pt < Pt max > where Pt max = (°- 059 + E^/Qll) GeV 2 . 
This cut, which is based on Monte Carlo signal studies, 
was designed to retain 90% of the signal. The efficiency 
after all cuts for IMD events was 79.6%. The final ef- 
ficiency for LNV events was 89.8% for V— A couplings 
and 79.9% for scalar couplings. Figures || and || show the 
pi distributions for right-sign and wrong-sign events re- 
spectively. Right-sign events are neutrino (anti-neutrino) 
events with an outgoing pT (/i + ); wrong-sign events are 
the opposite: neutrino (anti-neutrino) events with a out- 
going fi+ (fj,~). 

For right-sign events, the primary backgrounds that 
enter the IMD sample come from low hadronic energy 
neutrino charged-current interactions in the detector. 
These include quasi-elastic events, resonance events, and 
some small fraction of deep-inelastic scattering (PIS) 
events with a very small momentum transfer ^,0. Un- 
derstanding the background levels is essential for the 
IMD measurement, since the signal to background ratio 
for these events is one to eight. To accomplish this, we 
perform a full Monte Carlo simulation of low hadronic en- 
ergy neutrino processes. To simulate neutrino resonance 
production, we used a low-Q 2 higher-twist approxima- 
tion |i"l}| . We found this method more accurate in av- 
eraging over all low-multiplicity states than the single- 
pion production model from Rein and Sehgal fl2(] . Nu- 
clear effects such as Fermi motion [[l3| and Pauli sup- 
pression ||,[l4j were also applied to the Monte Carlo sim- 
ulations. The Monte Carlo was absolutely normalized to 
data DIS events with hadronic energies above 30 GeV 
for each running mode. The normalization sample con- 
tained 0.83(0.25) million neutrino (anti-neutrino) inter- 
actions with a mean energy of 140(120) GeV. 

The dominant systematic uncertainties for the right- 
sign events are related to the modeling of these low 
hadronic energy processes. Systematic errors include ef- 



TABLE I. Errors on IMD expected signal and LNV ex- 
pected background. Total statistical and systematic errors 
reflect errors from full parameter fit, which take into account 
correlations between errors. 



Category 


IMD (%) 


LNV (%) 


Statistical Error 


±6.7 


±13.0 


Muon Energy Scale 


±1.0 


±2.2 


Hadron Energy Scale 


±0.3 


±0.7 


Angle Smearing 


±0.6 


±1.4 


Normalization 


±6.4 


±2.6 


Quasi-Elastic Cross-Section 


±1.0 


±0.6 


Pauli Suppression 


±8.5 


±2.2 


Beam Impurities 


N/A 


±0.7 


Charge Identity 


<0.1 


±1.5 


Radiative Corrections 


±1.0 


±1.0 


Total Systematics (fit) 


±8.2 


±4.4 



fects from muon energy and angular resolution, back- 
ground cross-section uncertainties, Pauli suppression, 
and MC normalization. In addition, we take into account 
radiative correction errors which affect the IMD cross- 
section. A complete list of systematic errors is shown in 
Table | 

The validity of the background modeling was checked 
directly against the data by looking at the right-sign, 
anti- neutrino process v ^ + N — > fi + + N'. This partic- 
ular configuration selects only background events, and 
thus is an ideal platform to test the data to Monte Carlo 
agreement and systematics. A fit to the anti-neutrino pi 
distribution (Fig. ^ is performed where the backgrounds 
are allowed to vary within the uncertainties shown in 
the first column of Table [| The fit gives an excellent 
X 2 /d.o.f. of 44.9/50, indicating that the background es- 
timate agrees well with the anti-neutrino data within the 
systematic uncertainties. 

Having verified the size and spectrum of the back- 
ground, a fit to the neutrino data is performed to extract 
the IMD signal. The fit includes the previously men- 
tioned backgrounds plus an IMD signal contribution with 
the proper pi distribution. As before, the backgrounds 
are allowed to vary within the uncertainties shown in the 
first column of Table Q. As shown in Fig. [|, the data are 
well described by the combination of an IMD signal at 
low pi plus the background. From the fit, we extract a 
total of 1050 ± 139 IMD events, where 1311 events were 
expected based on Standard Model predictions, taking 
into account radiative cor rections |l| (see Table ||). 

The differential cross-section for IMD can be written 
as: 

^ = a -E u -(l-r) (3) 
ay 

where y = E^/Ey, r = m^/s, a = 2 "^ Gf , and s 
is the center-of-mass energy of the interaction. For in- 



TABLE II. Signal extraction from Monte Carlo background. 



Type 


v Mode / /i Charge 


Data 


Fit Results 


IMD Signal 


V 1 lT 


11792 


1050 ± 139 


LNV Signal (V-A) 


v 1 \T 


24 


0.6 ±3.3 


LNV Signal (scalar) 


V 1 fl" 


24 


-0.6 ±3.3 
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FIG. 3. Transverse momentum distributions for data 
(crosses) and Monte Carlo (solid) for neutrino events (left) 
and anti-neutrino events (right) with a wrong-sign muon. 
A LNV signal would appear as an excess of events in 
anti-neutrino mode. 

verse muon decay, the NuTeV measurement for the IMD 
asymptotic cross-section (E v ^> m^) is: 

(To = (13.8 ± 1.2 ± 1.4) x 10~ 42 cm 2 /GeV (4) 

where the first error is statistical, and the second is sys- 
tematic. The average neutrino energy for the IMD events 
sampled in the NuTeV experiment is 130 GeV. This mea- 
surement is in agreement with the theoretically predicted 
value of 17.2 x 10~ 42 cm 2 /GeV and is also consistent 
with the CHARM II measurement of (16.5 ±0.9) x 10~ 42 
cm 2 /GeV @. 

By requiring that the muon charge not match the neu- 
trino running mode (wrong-sign events), the analysis im- 
mediately becomes sensitive to lepton number violation. 
The dominant backgrounds in this case arise from beam 
impurities and muon charge mis-identification. Beam im- 
purities come mainly from charmed meson decays and 
decays of wrong-sign hadrons produced in secondary in- 
teractions p7| . Beam impurities constitute about 72% 
of the total LNV background. Charge mis-identification 
backgrounds are often associated with (5-ray production 
or multiple scattering of the muon in the toroid spectrom- 
eter. These backgrounds can be greatly reduced by im- 
posing quality cuts on the muon track in the toroid spec- 
trometer. The total fraction of charge mis-identification 
is 0.06% for anti-neutrino running mode. This source 
contributes 14% of the LNV background. Finally, there 
exists an irreducible background from V e e~ — ► y^V^ 
which contributes about 14% of the LNV background. 



The generic expression for the differential LNV cross- 
section is given by: 



da 
dy 



X^£l[A v .y(y-r)+A s -(l-r)] (5) 



where A represents the strength of the interaction and 
Ay and As determine whether the reaction is V— A or 
scalar. Integrating over all allowed values of y, and nor- 
malizing to the Standard Model IMD cross-section, al- 
lows the LNV cross-section to be written as: 



oiy^e -> fi v e ) 



A ■ [A v ■ ( 



l+r/2 x 



A, 
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We can make a consistency check on the background 
estimation by looking at the neutrino process + N — > 
/./+ + N'. Momentum distributions of these events are 
shown in Fig. ||. A total of 28 data events were seen in 
this sample where 23.5±3.7 (stat. + sys.) were expected, 
consistent with the background estimate. 



Looking in the LNV signal channel 



yields a total of 24 data events. A fit of the pi distri- 
bution to only background sources yields an acceptable 
X 2 /d.o.f. of 2.5/5, showing no indication of a LNV sig- 
nal. Including a possible LNV signal in the fit yields a 
total LNV contribution of 0.6±3.1±1.1 events for a V— A 
coupling and — 0.6±3.1±1.1 events for a scalar coupling. 
As shown in Table |[ the LNV analysis is dominated by 
statistical uncertainty. These fit results can be recast in 
the form of 90% C.L. limits on the LNV cross-section as 
a function of {Ay — As) / (Ay + As), as shown in Figure ||. 
If we assume a pure V— A coupling, this yields A < 1.7% 
while a scalar coupling yields a limit of A < 0.6%. 

This limit is currently the most stringent limit ob- 
tained directly from neutrino-electron scattering. Pre- 
vious results limited the pure V— A coupling to below 
5% jfj| . The LAMPF experiment has set an earlier limit 
from muon decay rates at < 1.2% |l!| for pure V— A 
couplings. 

In conclusion, NuTeV has performed a measurement of 
the inverse muon decay cross-section and a direct search 
for lepton number violation. The IMD asymptotic cross- 
section is measured to be (13.8 ± 1.2 ± 1.4) x 10~ 42 
cm 2 /GeV. The LNV search limits the strength of the 
interaction to be < 1.7% for V— A and < 0.6% for scalar 
couplings. 

This research was supported by the U.S. Department of 
Energy and the National Science Foundation. We thank 
the staff of FNAL for their contributions to the construc- 
tion and support of this experiment during the 1996-1997 
fixed target run. 
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FIG. 4. Limit on the lepton number violation process 
Vfj,e~ — ► /i~77 e as a function of the scalar and vector couplings 
(Ay — As)/(Av + As)- The non-shaded region is excluded 
by this result. 
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